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ABSTRACT 



^ ■ Recent progress in the theory of the electroweak phase transition is 

i/-) . discussed. It is shown, that for the Higgs boson mass smaller than the masses of 

W and Z bosons, the phase transition is of the first order. However, its strength is 
approximately 2/3 times less than what follows from the one-loop approximation. 
This rules out baryogenesis in the minimal version of the electroweak theory with 
light Higgs bosons. The possibility of the strongly first order phase transition in 

Q_i! the theory with superheavy Higgs bosons is considered. 

Ph. We show that if the Yang-Mills field at high temperature acquires a magnetic 

mass ~ g 2 T, then the infrared problem and the problem of symmetry behavior 
at high temperature effectively decouple from each other, no linear terms appear 
in the effective potential in all orders of perturbation theory and the symmetry 
in gauge theories at high temperatures is actually restored. Even though the 
last statement was never questioned by most of the authors, it was extremely 
difficult to come to a reliable conclusion about it due to the infrared problem in 
thermodynamics of non-Abelian gauge fields. 

The phase transition occurs due to production and expansion of critical 
bubbles. A general analytic expression for the probability of the bubble formation 
is obtained, which may be used for study of tunneling in a wide class of theories. 
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1 Introduction 



This talk is based on the results obtained in our papers with Michael Dine, Patrick Huet, 
Robert Leigh and Dmitri Linde |lj]. It contains also some more recent results on the in- 
frared problem in the electroweak theory, on the existence of symmetry restoration at a high 
temperature and on the possibility of the first order phase transition in the theory with 
superheavy Higgs fields. 

The existence of the phase transition in the electroweak theories was discovered by David 
Kirzhnits twenty years ago 0. A detailed theory of the phase transition was proposed in 
1974 by three groups of authors independently (by Weinberg, Dolan and Jackiw and by 
Kirzhnits and Linde ||), and soon the theory of the electroweak phase transition became 
one of the well established ingredients of modern cosmology. Surprisingly enough, the theory 
of this phase transition is still incomplete. 

In the first papers on this problem it was assumed that the phase transition is of the 
second order |2|, |3|]. Later Kirzhnits and Linde showed Q that in the gauge theories with many 
particles, and especially with particles which are much more heavy than the Higgs boson 
<p, one should take into account corrections to the high temperature approximation used in 
[0. |J. These corrections lead to the occurrence of cubic terms ~ g 3 (p 3 T in the expression 
for the effective potential V(<p, T). As a result, at some temperature, V acquires an extra 
minimum, and the phase transition is first order [|J. Such phase transitions occur through 
the formation and subsequent expansion of bubbles of the scalar field <fr inside the symmetric 
phase <p = 0. A further investigation of this question has shown that the phase transitions 
in grand unified theories are always strongly first order ||. This realization, as well as the 
mechanism of reheating of the universe during the decay of the supercooled vacuum state 
suggested in ||, ||, played an important role in the development of the first versions of 
the inflationary universe scenario [7[. (For a review of the theory of phase transitions and 
inflationary cosmology see Ref. ||.) 

For a long time it did not seem likely that the electroweak phase transition could have 
any dramatic consequences. Even though the possibility of a strong baryon number violation 
during the electroweak phase transition was pointed out fifteen years ago 0, Kj, only after 



the paper by Kuzmin, Rubakov and Shaposhnikov [I Ij was it realized that such processes 
do actually occur and may erase all previously generated baryon asymmetry of the universe. 

Recently, the possibility that electroweak interactions may not only erase but also produce 
the cosmic baryon asymmetry has led to renewed interest in the electroweak phase transition. 



A number of scenarios have been proposed for generating the asymmetry |12| - |L9| . All of 
them require that the phase transition should be strongly first order since otherwise the 
baryon asymmetry generated during the phase transition subsequently disappears. In all of 
these scenarios the asymmetry is produced near the walls of the bubbles of the scalar field <fi. 
Therefore it is necessary to make a much more thorough analysis of the electroweak phase 



transition than the analysis which is necessary for an approximate calculation of the critical 
temperature. 

We will say that the phase transition is strongly first order if the ratio of the Higgs field 
inside the bubble to the temperature T is larger than one, since otherwise the baryon 
asymmetry will be washed out by nonperturbative effects. 

This condition was used in [12|, ^] to impose a strong constraint on the Higgs mass in 
the minimal version of the electroweak theory, m# < 42 GeV. This, of course, already 



contradicts the present experimental limits m# > 57 GeV pl| . However, more careful 
consideration of various theoretical uncertainties indicated that the constraint might be 
somewhat weaker, permitting irtu up to 55 GeV, or possibly higher 2^j . In multi-Higgs 
models PUL IT5J, the limits are substantially weaker. 



Before one can discuss details of the process of baryogenesis, it is necessary to check 
that the results of our investigation of the phase transition are reliable. This is not a trivial 
issue even in the minimal electroweak theory. Indeed, as stressed in Refs. 0, ||, each new 
order of perturbation theory at finite temperature may bring a new factor of g 2 T ' jvn ~ gT/(p 
for the theories with gauge boson masses m ~ g<p. This means that the results of the one 
loop calculations may become unreliable at < gT. Thus, it became very desirable to go 
beyond the one-loop approximation. 

An example of such an approach is given by the self-consistent approximation elaborated 
in 0] . In this approximation, instead of the mass of a particle at zero temperature one uses 
its temperature-dependent mass, taking into account the contribution from the polarization 
operator. This method made it possible, in particular, to overcome unphysical difficulties 
related to imaginary masses of scalar particles at small 0. 

Recently this approach was reinvented by many authors. Some of the recent results 
obtained by this method were quite surprising. For example, it was claimed that higher 
order corrections lead to the appearance of a term in the effective potential ~ —g 3 <f)T 3 
[[23], |23j . This term is linear in 0; it is very large at small 0. Depending on its sign, it either 
may remove the local minimum of V(4>, T) created by the cubic term ~ — g 3 3 T, or it may 
make this minimum much more deep. 

Our investigation of this problem shows that if one is careful with counting of Feynman 
diagrams, neither positive nor negative linear terms ~ g 3 <j)T 3 appear in the effective potential 
[pi] . Even though now the authors of Refs. G3, EH] agree that the linear terms ~ g 3 <j)T 3 are 



absent, we will repeat our main arguments here, since these arguments may allow us to do 
much more than just say that there are no linear terms in order g 3 . A generalization of these 
arguments allows us to formulate the conditions under which one can show, despite some 
uncertainties with higher order corrections, that the expectation value of the scalar field 
at high temperatures actually disappears, = 0. Note, that this would be impossible in the 
presence of linear terms of any magnitude and sign. 



However, higher order corrections do lead to a significant modification of the one-loop 
results. They lead to a decrease of the cubic term g 3 (f) 3 T by a factor 2/3 JIJ.Q This effect 
decreases the ratio (p/T at the point of the phase transition by approximately the same factor 
2/3. This makes baryogenesis virtually impossible in the context of the minimal standard 
model with m w > m H > 57 GeV. 

All results discussed above were obtained in the context of the theories with small coupling 
constants and light Higgs fields. However, we do not really know whether the Higgs boson is 
light or very heavy. If the Higgs boson is superheavy, irtn ^ 10 3 GeV, the phase transition 
may become strongly first order. Even though this possibility is extremely speculative, it 
may lead to important consequences. Therefore we will discuss it in this paper. 

Assuming that one knows the shape of the effective potential at small 0, one should still 
work hard to determine the ratio <p/T at the point of the phase transition. One needs to 
know at what temperature the transition actually occurs, and some details of how it occurs. 
At very high temperatures the effective potential of the Higgs field, V(<f), T), has a unique 
minimum at the symmetric point = 0. As the temperature is lowered, a second minimum 
appears. At a critical value T c , this second minimum becomes degenerate with the first 
one. However, the phase transition actually occurs at a somewhat lower temperature, due 
to the formation of bubbles of true vacuum which grow and fill the universe. The usual 
way to study bubble formation is to use the euclidean approach to tunneling at a finite 
temperature [27]. One should find high-temperature solutions, which describe the so-called 



critical bubbles. Then one should calculate their action, which leads to an exponential 
suppression of the probability of bubble formation. Typically, these calculations are rather 
complicated, and analytic results can only be obtained in a few cases. One of these is the 
thin wall approximation, which is valid (as in the case of transitions at zero temperature) 
if the difference in depth of the two minima of V(4>, T) is much smaller than the height 
of the barrier between them. In this case the radius of the bubble at the moment of its 
formation is much larger than the size of the bubble wall, and the properties of the bubble 
can be obtained very easily. However, the thin wall approximation in our case leads to an 
underestimate of the tunneling action by a factor of two. Fortunately, we were able to obtain 
a simple analytic expression which gives the value of the euclidean action for theories with 
an effective potentials of a rather general type, V(<p, T) = a<j) 2 — 60 3 + c0 4 . We hope that this 
result will be useful for a future investigation of bubble formation in a wide class of gauge 
theories with spontaneous symmetry breaking. 

On the other hand, validity of the standard assumption that the phase transition occurs 



due to formation of critical bubbles should be verified as well. Kolb and Gleiser 281 and 



more recently, Tetradis |2I| have argued that the phase transition may occur by a different 
mechanism, the formation of small (subcritical) bubbles. If this is the case, the transition 

2 A similar result was obtained also by Carrington |]26| . However, her original results were different from 
ours approximately by a factor of two. Consequently, they lead to an impression that modification of the 
one-loop results should lead to an increase of the strength of the first order phase transition. At present, 
there is no disagreement between our results. 



is completed earlier and by a different mechanism than in the conventional picture. While 
this idea is very interesting, we will argue (see also BO] and the talk of Anderson at this 



Conference) that it is only relevant in cases where the transition is very weakly first order 
and the euclidean action corresponding to critical bubbles is not much larger than one. This 
is not the case for the strongly first order phase transitions, where the relevant value of the 
euclidean action at the moment of the transition is S ~ 130 — 140. 



2 The Phase Transition 

Let us consider the form of the effective potential at finite temperature. Contributions of 
particles of a mass m to V((j),T) are proportional to m 2 T 2 , m?T and m 4 ln(m/T). We will 
assume that the Higgs boson mass is smaller than the masses of W and Z bosons and the 
top quark, m# < mw,mz,m t . Therefore we will neglect the Higgs boson contribution to 
V(cf>,T). 

The zero temperature potential, taking into account one-loop corrections, is given by || 

V = -^0 2 + V + 2Bv 2 J 2 - jW + 50 4 ln(^) . (1) 



Here 

B =64^ (2rr4 + m| _4m * 4) ' (2) 

v = 246 GeV is the value of the scalar field at the minimum of Vq, A = fi 2 /v 2 , m 2 H = 2/i 2 . 
Note that these relations between A, //, v and the Higgs boson mass m H , which are true at 
the classical level, are satisfied even with an account taken of the one-loop corrections. This 
is an advantage of the normalization conditions used in ||. An expression used in |Z3| is 
equivalent to this expression up to an obvious change of variables. 

At a finite temperature, one should add to this expression the term 

v T = ^(Qi-(yw) + 3i-(yz) - Qi + (y t )) , (3) 

where yi = Mi<fi/v T, and 

I T (y) = ± f°° dxx 2 \n(lTe~^ T ^) . (4) 

Jo 

The results of our work are based on numerical calculation of these integrals, without making 



any specific approximations [22]. However, in the large temperature limit it is sufficient to 



use an approximate expression for V((f),T) |fj, |23 



V&, T) = D(T 2 - T 2 )<P 2 - ETf + ^0 4 . (5) 



Here 



D = -^(2m 2 w + m 2 z + 2m 2 t ), (6) 

E = -^(2m 3 w + m|) ~ 1(T 2 , (7) 



T 2 = ^ 2 - 4B, 2 ) = ^(m| - 8Bvl) 






X T = X 3 —( 2 mt v ln^ + mi,ln^-4mtln^-) , (9) 

16vr 2 ^ ^ w a B T 2 z a B T 2 ' a F T 2 J ' v ; 

where lna B = 21n47r — 27 ~ 3.91, lna^ = 21ri7r — 27 ~ 1.14. 

It will be useful for our future discussion to identify several 'critical points' in the evolution 
of\/(0,T). 

At very high temperatures the only minimum of V((f>, T) is at — 0. A second minimum 
appears at T — T 1; where 

^ = l-9E>°/8\ Tl D ' (10) 

The value of the field in this minimum at T — 7\ is equal to 

3£T\ 



2A 



:i^ 



ii 



The values of V(0, T) in the two minima become equal to each other at the temperature T c , 
where 

T < = T^o ■ < 12 > 

At that moment the field in the second minimum becomes equal to 

■ 2EZ (13) 



The minimum of V((f>, T) at = disappears at the temperature T , when the field in the 
second minimum becomes equal to 

3ET ° (14) 



A 



To 



3 Infrared Problems and Reliability of the Perturba- 
tion Expansion 

In our previous discussion, we have considered only the one loop corrections to the effective 
potential. In this section we discuss the role of higher order corrections. 
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It is well known that, in field theories of massless particles, perturbation theory at finite 
temperature is subject to severe infrared divergence problems. For small values of the scalar 
field, the gauge bosons (and near the phase transition, the Higgs boson) are nearly massless; 
as a result, as was pointed out in the early work on this subject j|, U, one cannot reliably 
compute the effective potential for very small <fi. The problem is that the higher order 

/ a 2 T\ N 

corrections in coupling constants may contain terms of the type of (^-) • As a result, 
higher order corrections go out of control for m < g 2 T. For scalar particles this happens 
near the critical point only. Indeed, scalar particles have masses m ~ gT ^> g 2 T in the high 
temperature limit. However, if one takes into account gauge invariance, it can be shown that 
"magnetic" components of vector particles cannot acquire any contribution to their "masses" 
larger than g 2 T. 

At this point one should be more precise. The Green function of the vector field is 
singular at k 2 ~ g 2 T 2 . In this sense one may speak about the vector field mass ~ gT. 
However, the Green function of the vector field at a finite temperature does not have a 
simple pole singularity. For example, in addition to the singularity at k 2 ~ g 2 T 2 , the Green 
function of a photon has a singularity at k = 0, k — > 0. It is this singularity that is 
responsible for all infrared problems in quantum statistics of gauge fields, since the Green 
functions at k$ = give the leading infrared divergent contribution to thermodynamical 
sums H [33|, [34) . Investigation of the infrared problem in gauge theories without spontaneous 
symmetry breaking has shown that the "magnetic mass" (corresponding to the limit k = 
0, k — ► 0) may appear in the non-Abelian theories, but it cannot be larger than 0(g 2 )T. 

Thus, in the absence of spontaneous symmetry breaking, or at < gT, when the magnetic 
mass of the vector particles become smaller than g 2 T, perturbative results may become 
unreliable. One may wonder, therefore, is it possible that some unusual contribution to the 
effective potential at <fi < gT may alter our results. 

Recently, in a very interesting paper, Brahm and Hsu found that at small <ft, higher order 
corrections to the scalar field contribution to the effective potential may produce a large 
negative linear term — g 3 0T 3 , which eliminated any trace of a first order transition. 

On the other hand, Shaposhnikov considered higher order corrections to the vector par- 
ticle contribution to V((f>,T) and found a large positive term + g 3 0T 3 which made the phase 
transition strongly first order (4>/T > 1) even for m# ~ 64 GeV 



We will show that neither positive nor negative linear terms appear in the expression 
for V((j),T) if one studies higher order corrections paying particular attention to the correct 
counting of Feynman diagrams. [||] . (Additional information on this problem is contained in 
the talks by Michael Dine and Robert Leigh at this conference.) 

We will consider here for simplicity the contribution of the scalar particles and the W 
bosons only; adding the contribution of Z bosons is trivial. As we have already noted, for 
questions of infrared behavior, fermions may be ignored. Coulomb gauge, V • W = 0, is 



particularly convenient for the analysis, though the problem can be analyzed in other gauges 
as well. In this gauge, the vector field propagator D^ u after symmetry breaking (and after 
a proper diagonalization) splits into two pieces, the Coulomb piece, D 00 , and the transverse 
piece, Dij. For non-zero values of the discrete frequency, u n = 2nnT, the Coulomb piece 
mixes with the 'Goldstone' boson. However, for the infrared problems which concern us here, 
we are only interested in the propagators at zero frequency. For these there is no mixing. 
One has [§ 

D 0O (u = 0,k) = gr - (15) 

k 2 + m w {<p) 
and 

A> = o, k) = r 1 —— Pij (k) , (16) 

k 2 + m w {<p) 

where Py = 5ij y-. The mass of the vector field W at the classical level is given by 

m w — gVo/Z- Propagators of the Higgs field (f) and of the 'Goldstone' field \ in this gauge 
are given by 

D^) = j^— 2 , (17) 

k 2 + m^ 

D x (k) = i • (18) 

k 2 

Let us review several ways of obtaining the standard one-loop expression for the cubic 
term in the effective potential, eq. (|5|). The most straightforward is to carefully expand eq. 
(H) for the effective potential in yw = -^fr = fr- Indeed, the contribution of VF-bosons to 
the effective potential at T > mw(4>) is given by 

wvr ' ; V 90 24 12tt 

= 2x3xf-^ + ^f-^T + ...). (19) 

^ 90 96 96tt > y ' 

Here the expression in brackets coincides with the contribution of a scalar field with mass mw] 
the factor 2 appears since there are two jy-bosons with opposite charges, while the factor 3, 
which will be particularly important in what follows, corresponds to the two transverse and 
one longitudinal degrees of freedom with mass mw- 

Alternatively, we can obtain the cubic term by looking directly at the one-loop Feynman 
diagrams. For this purpose, it is only necessary to examine the zero frequency contributions. 
Certain diagrams containing four external lines of the classical scalar field naively give a 
contribution proportional to g 4 4 ; the cubic term arises because the zero frequency integrals 
diverge for small mass as T /mw ~ T/g<j). 

Consider, in particular, the zero frequency part of the expression for the one loop free 
energy in momentum space. It is simplest to compute the tadpole diagrams for dV/dcj) and 
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afterwards integrate with respect to 0. The transverse gauge bosons give a contribution 
dV tI g 2 <j)T f d 3 k 1 g 2 (j) T 



«, 2x f !wfw^r- 2x ^ ^' (20) 

where, by keeping only the zero frequency mode, we have dropped terms which are analytic 
in m 2 . The Coulomb lines give half the result of eq. fl2"0"|). Integration of the total vector 



field contribution correctly represents the cubic term in (19). 



A complete gauge boson contribution to the tadpole, including the non-zero frequency 
modes, is [§J 

^f T) = 2x3x^(f-^ + ...) = 2x3x^(T 2 -^ + ...). (21) 



48 V 7i ' 48 ^ 2tt 

One can easily check that integration of this expression with respect to gives eq. (|H]) . 

With these techniques, we are in a good position to study higher order corrections to 
the potential. The authors of Refs. [[24, ^] found a linear contribution to the potential by 



substituting the mass found at one loop back into the one loop calculation. The effective 
masses-squared of both scalar particles and of the Coulomb field contain terms of the form 
~ g 3 T0, which, upon substitution in (34), give linear terms. But this procedure is not always 
correct. It is well known that the sum of the geometric progression, which appears after 
the insertion of an arbitrary number of polarization operators 11(0, T) into the propagator 
(k 2 + m 2 )" 1 , simply gives (k 2 + m 2 + 11(0, T)) _1 . Therefore one can actually use propagators 
(k 2 + m 2 + n(T)) -1 , which contain the effective mass-squared m 2 + 11(0, T) instead of m 2 . 
However, this trick with the geometric progression does not work for the closed loop diagram 
for the effective potential, which contains \n(k 2 + m 2 ). A naive substitution of the effective 
mass squared m 2 + 11(0, T) instead of m 2 into ln(fc 2 + m 2 ) corresponds to a wrong counting 
of higher order corrections. 

This does not mean that there is no regular way to make this trick for the effective 
potential. One may add and subtract from the Lagrangian the term — ~ 2 I1(0(T),T) to 
the Lagrangian, and a similar term for the vector field as well. Here n(0(T),T) is the 
polarization operator with an account taken of all daisy and superdaisy diagrams, 4>{T) is 
a classical field, not an operator. Then the effective mass (at zero momentum) becomes 
renormalized, m 2 — > m 2 + n(0, T), but one should add some extra diagrams containing the 
insertion — | 2 n(0(T),T). These diagrams were not considered in f23| , P5fl . Here one can 
clearly understand the difference between calculating effective potential and tadpoles. If 
one uses this method to calculate the tadpole diagrams, insertions of — \ 2 n(0(T),T) self- 
consistently cancell the diagrams which would generate extraneous polarization operator 
corrections to already corrected effective mass m 2 (0(T),T) = m 2 + n(0,T). This gives us 
the standard prescription of simply substituting m 2 (<j)(T),T) = m 2 + n(0,T) instead of m 2 
in all propagators (k 2 + m 2 )~ 1 . 



Thus, the simplest way to take into account high temperature corrections to masses of 
vector and scalar particles without any problems with combinatorics is to compute tadpole 
diagrams for ^; these are then trivially integrated to give the potential. One can easily check 
by this method that no linear terms appear in the expression for V(<p, T). Indeed, at a given 
temperature and effective mass, the tadpoles are linear in (see e.g. equation ([21]) ) . To take 
into account the mass renormalization in the tadpoles, one should substitute the effective 
mass squared m 2 + 11(0, T) into the one-loop expression for the tadpole contribution; as we 
explained above (see also ([|J)), this is a correct and unambiguous procedure for tadpoles. 

dVt L = 2x 9^r r_fk 1 m) 

# 2 7 (2vr) 3 fr + m 2 + tt((f),k,T) ' 

This expression could lead to a linear term in the effective potential only if the integral 
would behave as _1 in at small — > 0. However, in our case this is impossible, since the 
polarization operator calculated in [F23], |25| at small k and is nonnegative, and the integral 
converges to a constant. Therefore its subsequent integration with respect to 0, which gives 
the correction to the effective potential, is quadratic in 0, i. e. it does not contain any linear 
terms. 

We should emphasize that our approach effectively takes into account all diagrams con- 
sidered in |[24| , ||25|| , but with correct combinatorics. Note also, that this approximation works 



even if the polarization operator depends on the classical field 0.0 An apparent triviality of 
the investigation with the help of tadpoles is not due to its incompleteness, but due to the 
power and simplicity of this method, which was elaborated in |4j] . In particular, it was even 
unnecessary for us to use a particular expression for the polarization operator, as far as it is 
nonnegative. We will return to this issue shortly. 

Even though there are no linear terms ~ g 3 0T 3 , higher order corrections do have a 
dramatic effect on the phase transition. This effect is a modification of the cubic term. 

As we have shown above, the cubic term appears due to the contribution of zero modes, 
uj n = 2imT = 0. This makes it particularly easy to study its modification by high order 
effects. Indeed, it is well known that the Coulomb field at zero frequency acquires the Debye 
'mass', m 2 D = Hoo^n = 0, k — * 0) ~ g 2 T 2 . This leads to an important modification of the 
Coulomb propagator (|T5|): 

Ax)(fc-+0)= 2 \ (23) 

k z + m z D + mf v [(p) 

For the values of of interest to us, m 2 D 3> m 2 v ((f)). Thus, repeating the calculation of 
the cubic term, the Coulomb contribution disappears. However, the transverse contribution, 

3 We disagree with the recent claim made in [ p5[ that our approximation works only if the polarization 
operator II does not depend on </>, and that we do not take into account 'subleading' diagrams which give 
rise to dangerous linear terms obtained in |24|], p5[. Indeed, in Q we did not take into account diagrams of 
the type shown in Fig. 2; see their discussion below. However, these diagrams were not considered in pij, 
J25fl as well, since they do not lead to the linear terms in the order g 3 . 
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which is two times larger than the Coulomb one, is unaffected at this order, due to the 
vanishing of the 'magnetic mass' [[], |33|, |3"1| . As a result, the cubic term does not disappear, 
but it is diminished by a factor^ of 2/3: 

E = ^{2ml + m%). (24) 

This small correction proves to be very significant. Indeed, eqs. (|HJ), ( |T^D show that 
the ratio of the scalar field to the temperature at the moment of the phase transition is 
proportional to E, i.e. to the cubic term. Actually, the dependence is even slightly stronger, 
since for smaller E the tunneling occurs earlier. Even before the reduction of the cubic 
term was taken into account, the ratio <f)/T for rrin ^ 57 GeV was slightly less than the 
critical value <p/T « 1. The decrease of this quantity by a factor of 2/3 makes it absolutely 
impossible to preserve the baryon asymmetry generated during the phase transition in the 
minimal model of electroweak interactions with m# > 57 GeV. 

Are these results completely reliable? The effective coupling constant of interactions 
between W bosons and Higgs particles is g/2. In this case, a general investigation of the 
infrared problem in the non-Abelian gauge theories at a finite temperature suggests that the 
results which we obtained are reliable for > | T ~ T/3 ||, [^, 0. Thus, a more 
detailed investigation is needed to study behavior of the theories with m H > 10 2 GeV 
near the critical temperature, since the scalar field, which appears at the moment of the 
phase transition in these theories, is very small. However, we expect that our results are 
reliable for strongly first order phase transitions with > T/3, which is quite sufficient to 
study (or to rule out) baryogenesis in the electroweak theory. Recent investigation of higher 



order corrections to this theory indicates |36[ that these results may be reliable even down 
to (j) ~ gT/10. 

Finally, we would like to address a fundamental question: since the theory for <fi <C 
gT is infrared divergent, can we definitely establish that the symmetry is restored at high 
temperature, or is it possible that <ft always has some small, non-zero value? Indeed, if our 
approximation breaks down at <ft < gT, how do we know that the symmetry restoration 
actually takes place, i.e. <fi = at T > T ? 

To address this question, we can work far away from the critical point, at T — T ^> T . 
The best approach, as before, is to study all possible higher order tadpole diagrams, see Figs. 
1,2. 

There are two different classes of diagrams to be considered. External line of the scalar 
field may split either into two lines of the vector field, Fig. 1, or into two lines of vector 
field and one line of scalar field, Fig. 2. All diagrams of the first type can be represented 
as the trivial one-loop diagram, plus the diagrams with an arbitrary number of polarization 



4 There was also a claim that the cubic term disappears completely J37|, but recently this claim was 
withdrawn. 
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operator insertions. The simplest diagram of that type is shown in Fig. 1. The black circle 
stands for the exact polarization operator U((f),k,T), including all higher order corrections. 
The sum of all these diagrams gives us the one-loop diagram with an exact Green function 
of the vector field instead of the free field propagator. In other words, as usual, one must 
add polarization operator to the mass squared of the vector field. 

The behavior of 11(0, k, T) at k > g 2 T is known perturbatively. It leads only to high-order 
corrections to V((j),T), which do not contain any nonanalytic terms, such as the dangerous 
linear terms in 0. The only possible source of problems is our absence of knowledge of 
n(0, k,T) at < gT, k$ = 0, \k\ < g 2 T. Indeed, in this domain all higher order corrections 
to II((/>, k,T) are equally important. However, the consequences of this uncertainty may 
appear not very significant. 



Indeed, consider again the most dangerous part of the tadpole diagram, eq. (|20|), and 

add U((p,k,T) to m 2 w ((p): 

dV^ = 2x g^T r d*k 1 

d<p 2 J (2vr)3 fc2 + ^! + n(0,A;,T) ' l ' 

The part of the integral in the domain of uncertainty, \k\< g 2 T, is given by 

g 2 <j)T [9 2 t k 2 dk , . 

1 (26) 



2vr 2 Jo k 2 + sff + U((f),k,T) ' 

On dimensional grounds, one expects that at g<p <C g 2 T and \k\ <C g 2 T the polarization 
operator has some value of the order g A T 2 || |33], [34] . This follows from the fact that the 
most infrared divergent part of the theory corresponds to the three-dimensional theory, with 
g 2 T being the only mass (or coupling constant) scale. 

One cannot exclude a possibility that the polarization operator is proportional to — g 2 \k\T, 
or it is of the order g 4 T 2 , but is negative, and the integral in ( |2"6] ) diverges in the limit 
— > 0. One should keep this possibility in mind, since it may lead to interesting and 
unusual consequences, like Bose-condensation or even crystallization of the Yang-Mills field 
at high temperature |33|| . Indeed, the tadpole integral has a simple interpretation in terms 
of integration over the occupation numbers of bose fields. Large contribution to this integral 
may be interpreted as a result of Bose condensation of particles in a state with a nonvanishing 
momentum. In our case, this effect may also lead to absence of a complete symmetry 
restoration at T > T . Note, that this effect may occur only at < gT, as we anticipated. 
However, it is not quite clear that this possibility is physically viable. Whereas occupation 
numbers may be large, they can hardly be negative. 

The standard (and most conservative) assumption is that < n(0 = 0, k = 0) < g 4 T 2 
[^, [53|, |33J. This corresponds to generation of a magnetic mass < m 2 < g A T 2 . In such case 
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we may estimate the corresponding integral at small as 

g 2 <PT [9 2 t k 2 dk 



2ir 2 Jo k 2 + 0(l)g 4 T 



2 



g^T 2 . (27) 



This term after integration over does not give any linear terms in 0. It gives just a small 
correction to the quadratic part of the effective potential, AV ~ g A <f) 2 T 2 . Such corrections 
do not alter our conclusions concerning symmetry restoration at high temperature. Note, 
that this term corresponds to the sum of all most dangerous contributions to the tadpole 
diagrams of the type of Fig.l, to all orders in g 2 . 

Now let us consider the diagrams which contain internal lines of scalar field, Fig. 2. These 
diagrams may be dangerous near the critical point, where the scalar fields are almost massless, 
but they are much less dangerous than the diagrams of the first class at the temperature 
much higher than critical. The reason is that at high temperature the scalar field acquires a 
large mass m 2 ~ g 2 T 2 ^> g 4 T 2 . Therefore scalar particles by themselves do not lead to any 
infrared problems outside of a small vicinity of the critical point. The presence of such heavy 
particles effectively cuts infrared divergencies in the diagrams with vector particles as well. 
One can easily check that the diagram with one vector loop, Fig 2a, gives the contribution 
g 4 (p 2 T 2 to the effective potential, the diagram with two vector loops gives g 5 2 T 2 , the diagram 
with three vector loops gives g 6 <f) 2 T 2 . Starting with this diagrams, infrared problem becomes 
manifest in that each diagram of this type with higher number of vector loops gives the 
contribution of the same order g 6 <p 2 T 2 . 

Thus, the infrared problem in thermodynamics of gauge fields does not permit us to 
calculate the effective potential at high temperature to all orders of perturbation theory. 
The diagrams Fig. 2 contain uncertainties at the level of g 6 <p 2 T 2 ; the diagrams Fig. 1 
contain uncertainties at the level of g A <j) 2 T 2 . However, neither of these diagrams produce 
linear terms in 0, unless the Green functions of a massless Yang-Mills field has a pathological 
behavior at large temperature. As for the quadratic terms, they can be calculated at least 
with an accuracy up to g 3 <p 2 T 2 , or maybe even up to g 4 lng <p 2 T 2 . This is quite sufficient to 
calculate the critical temperature and to make a conclusion that at the temperature higher 
than critical the scalar field vanishes. 

These considerations indicate that the situation with the phase transitions in the non- 
Abelian gauge theories is probably the same as in the standard case: infrared problems may 
prevent a simple description of the phase transition in a small vicinity of the critical point 
(unless the phase transition is strongly first order), but everywhere outside this region, the 
symmetry behavior of gauge theories can be described in a reliable way. 
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4 First order phase transitions with superheavy Higgs? 

In our previous investigation we neglected the contribution of Higgs bosons to the one-loop 
effective potential. The reason was very simple: We have seen that the increase of the Higgs 
boson mass decreases the strength of the phase transition. This can be easily understood by 
considering a model of a single scalar field with the effective potential 



V = -£<P + -0 4 = -^0 2 + -0 4 (28) 



: I ( > -h — ( ' 

2 r 4 



Near the point of the phase transition, where the high-temperature approximation works well, 
one may investigate the symmetry behavior in this theory in a self-consistent approximation 
suggested in Q, where only cactus diagrams should be evaluated. In this approximation 
the effective mass of the scalar field and the first derivative of the effective potential in its 
extremum are simply related to each other: 

m 2 (T, 0) = 3A0 2 - /i 2 + n(T, m(T, 0)) , (29) 

and 

-^ = H(T,m(T,0)). (30) 

d(p 

Here V 1 is the one-loop contribution to the effective potential, 
ldV 1 3A r™ k 2 dk 



2vr27 ° yP + m^rWexp V^f^ 



(31) 



At the minimum of the effective potential with 4>(T) ^ 

dV l 



+ A0 d - /i 2 = , (32) 

which, together with Q2"5D and Q3"0D, gives m? = 2A0 2 (T) and 

1 dV 1 3A [°° k 2 dk XT 2 ( 3V2X<p{T) \ 



*** 2n2J ° y/F + 2A0^(T) (exp v^W _ i) 4 V ^ 

(33) 
On the other hand, the local minimum of V((f>) at <fi = disappears when m(T, 0) = 0. This 
gives the critical temperature 

T = 2v, (34) 

where v = fi/yX = 246 GeV, v > 4>(T). For small A the last term in eq. fl3"3"| ) is small at 
T m T , which implies that the phase transition is weakly first order. Actually, we cannot 



even say from eq. (|33D whether the phase transition is second order or weakly first order, 



since near the critical point the higher order corrections are large 

14 



Let us nevertheless use eq. (^) to make a bold estimate of conditions under which the 
phase transitions could be strongly first order. From eqs. fl3"5|), fl3"l| ) h follows that the jump 
of the scalar field at the critical temperature (assuming that T is close to the temperature 



of the phase transition) is given by "^ . It is larger than the temperature T (which is the 
condition for baryogenesis) if A > 9, or, equivalently, 

m H > 10 3 GeV . (35) 

For obvious reasons, this estimate should not be taken as a serious indication of existence 
of strongly first order phase transitions and baryogenesis with superheavy Higgs bosons. 
However, the stakes are high, and the possibility to have a strongly first order phase transition 
and baryogenesis in the strong coupling regime with superheavy Higgs bosons (technicolour?) 
should not be overlooked. 



5 Bubble Formation 

In the previous section we noted that the two minima of V(4>, T) become of the same depth 
at the temperature T c , eq. (|T2|). However, tunneling with formation of bubbles of the field 
corresponding to the second minimum starts somewhat later, and it goes sufficiently fast to 
fill the whole universe with the bubbles of the new phase only at some lower temperature T 
when the corresponding euclidean action suppressing the tunneling becomes less than 130 - 
140 |n], ^, PS"] . In |[Jj (see also [p^j ) we performed a numerical study of the probability of 



tunneling. Before reporting our results, we will remind the reader of some basic concepts of 
the theory of tunneling at a finite temperature. 

In the euclidean approach to tunneling (at zero temperature) ||32|| , the probability of 
bubble formation in quantum field theory is proportional to exp(— S4), where S4 is the four- 
dimensional Euclidean action corresponding to the tunneling trajectory. In other words, S4 
is the instanton action, where the instanton is the solution of the euclidean field equations de- 
scribing tunneling. A generalization of this method for tunneling at a very high temperature 



27] gives the probability of tunneling per unit time per unit volume 



P~A(T)-exp(-^). (36) 

Here A(T) is some subexponential factor roughly of order T 4 ; S3 is a three-dimensional 
instanton action. It has the same meaning (and value) as the fluctuation of the free energy 
F = V((f)(x),T) which is necessary for bubble formation. To find S3, one should first find an 
0(3)-symmetric solution, 0(r), of the equation 

^ + 1 ^ = v'U) (37) 

dr 2 r dr 
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with the boundary conditions <p(r = oo) = and d(j)/dr | r= o = 0. Here r = Jxf] the Xj are 
the euclidean coordinates, i = 1,2,3. Then one should calculate the corresponding action 




S 3 = ^J q r 2 dr[-\-^\ +V(<j>(r),T)}. (38) 

Usually it is impossible to find an exact solution of eq. (|37|) and to calculate S3 without 



the help of a computer. A few exceptions to this rule are given in Refs. ||, |27| . One of these 
exceptional cases is realized if the effective potential has two almost degenerate minima, such 
that the difference e between the values of V(4>, T) at these minima is much smaller than the 
energy barrier between them. In such a case the thickness of the bubble wall at the moment 
of its formation is much smaller than the radius of the bubble, and the action S3 can be 
calculated exactly as a function of the bubble radius r, the energy difference AV and the 
bubble wall surface energy Si: 

A-TT 

S 3 = ~— r 3 AV + A7rr 2 Si, (39) 

where 

Si = | d^2V{<t>,T) . (40) 

The radius of the critical bubble r can be found by finding an extremum of S^r). However, 
one must be very careful when using these results. Indeed, as can be easily checked, this 
extremum is not a minimum of the action, it is a maximum. Therefore, the action corre- 
sponding to the true solution of eq. fl3"T|) will be higher than the action of any approximate 
solution. As a result, one can strongly overestimate the tunneling probability by calculat- 
ing it outside the limit of validity of the thin wall approximation. In our case the thin wall 
approximation underestimates the tunneling action by a factor of two, i.e. it gives the proba- 
bility of tunneling about e~ 100 where the correct answer is e~ 200 . If the only thing one wishes 
to know is the time when the tunneling occurs, this error is not very important. It leads only 
to a few percent error in calculation of the temperature of the universe at the moment of the 
phase transition, since the tunneling action is extremely sensitive to even very small changes 
of the temperature. Thus, one may argue that the thin wall approximation is still useful. 
(See also the talk of Anderson at this Conference.) However, it is possible to determine the 
time of the phase transition with an accuracy of few percent without any study of tunneling: 
It is enough to say that the phase transition happens in the middle of the interval between 
T c and T . In order to obtain a complete description of the phase transition, including a 
correct shape of the bubble wall, one should go beyond the thin wall approximation. 

We would now like to obtain an analytic estimate of the probability of tunneling in the 
electroweak theory, which can be used for any particular numerical values of constants D, 
E and A^. As shown in Ref. p3|j , eq. (|[) in most interesting cases approximates V(<p, T) 
with an accuracy of a few percent. This by itself does not help very much if one must study 
tunneling anew for each new set of the constants. However, it proves possible to reduce this 
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study to the calculation of one function f(a), where a is some ratio of constants D, E and 
\t- In what follows we will calculate this function for a wide range of values of a. This will 
make it possible to investigate tunneling in the electroweak theory without any further use 
of computers. 

First of all, let us represent the effective Lagrangian L((f),T) near the point of the phase 
transition in the following form: 

L(cj>,T) = \{d^f - ^V + ET<f - ^ . (41) 

Here M 2 (T) = 2D(T 2 — T 2 ) is the effective mass squared of the field <fi near the point = 0, 
< is the value of the effective coupling constant At near the point of the phase transition 
(i.e. at T ~ T t , where T t is the temperature at the moment of tunneling). With a very good 
accuracy, the constants Aj, At 1; At c , Xt are equal to each other. 



Defining = ^f&, x = X/M, the effective Lagrangian can be written as: 

T,((b T\ = 

AE 2 T 2 



L (^ T ) = T^4^) 2 " \& + \& ~ f * 4 ] > ( 42 ) 



where 

AqM 2 

a = 2E*F 2 ■ (43) 

The overall factor 4j ^ y2 does not affect the Lagrange equation 

d 2 $ 2# ^ 3^ 2 1 . 

— - H = $ $ 2 + -a$ 3 . (44) 

dR 2 RdR 2 2 K J 

Solving this equation and integrating over d 3 X = M~ 3 d 3 x gives the following expression for 
the corresponding action: 

S 3 4.85 M 3 

T = ~^W X f{a) ■ (45) 

The function f(a) is equal |27j to 1 at a = 0, and blows up when a approaches 1. In the 



whole interval from to 1 this function, with an accuracy about 2%, is given by the following 
simple expression: 

"r 2.4 0.26 i , s 

/(«) = ! + - 1 + -— + - -j . (46) 



4 



1 — a 



a] 



In the vicinity of the critical temperature T , i.e. at AT = T~T <^ T , the action fl4o| ) 
can be written in the following form: 

S 3 38.8 D 3 / 2 /AT\ 3/2 /2A I?ATn 



T £ 2 



T) x/ (— ^T-)' (47) 



Using these results, one can easily get analytical expressions for the tunneling probability 
in a wide class of theories with spontaneous symmetry breaking, including GUTs and the 
minimal electroweak theory. 
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6 Subcritical Bubbles 



Despite our semi-optimistic conclusions concerning the infrared problem, it is still desirable 
to check that the whole picture of the behavior of the scalar field described above is (at least) 
self-consistent. This means that if the effective potential is actually given by eqs. (|), (|]), (H), 
(||), (p4|), then our subsequent description of the phase transition and the bubble formation 
is correct. Indeed, one would expect that the theory of bubble formation is reliable, since the 
corresponding action for tunneling S3/T is very large, S3/T ~ 130 — 140. However, recently 
even the validity of this basic assumption has been questioned. Gleiser and Kolb |28| and 



Tetradis |^J have argued that in many cases phase transitions occur not due to bubbles of 
a critical size, which we studied in section 3, but due to smaller, subcritical bubbles. We 
believe that these authors raise a real issue. However, we will now argue that this problem 
only arises if the phase transition is extremely weakly first order. 



The basic difference between the analysis of Ref. pq , |29| and the more conventional one 
is their assumption that at the time of the phase transition there is a comparable probability 
to find different parts of the universe in either of the two minima of V(<f),T). The main 
argument of Ref. J28], |2P| is that if the dispersion of thermal fluctuations of the scalar field 
< <ft 2 >~ T 2 is comparable with the distance between the two minima of V((j), T), then the 
field (f) "does not know" which minimum is true and which is false. Therefore it spends 
comparable time in each of them. According to [28fl, a kind of equilibrium between the 



domains of the two types is achieved due to subcritical bubbles with small action S 3 /T if 
many such bubbles may appear within a horizon of a radius if -1 . 

In order to investigate this question in a more detailed way, let us re-examine our own 
assumptions concerning the distribution of the scalar field prior to the moment at which the 
temperature drops down to T\, when the second minimum of V(<p, T) appears. According to 



IT]) , the value of the scalar field <p in the second minimum at the moment when it is formed 
is equal to 4>i — §3p"- For m# ~ 60 GeV (and taking into account the coefficient 2/3 in the 
cubic term) one obtains <fi\ ~ 0.4 T. Thermal fluctuations of the field <fi have the dispersion 
squared < <p 2 >= T 2 /12. (Note an important factor 1/12, which was absent in the estimate 
made in ||28| .) This gives dispersion of thermal fluctuations y/< 2 > ~ 0.3 T, which is not 



much smaller than 0i. 

However, as the authors of [pj| emphasized in their previous work |38| (see also p9|), the 



total dispersion < >~ T jYl is not an adequate quantity to consider since we are not 
really interested in infinitesimally small domains containing different values of fluctuating 
field 0. They argue that the proper measure of thermal fluctuations is the contribution to 
< cf) 2 > from fluctuations of the size of the correlation length £(T) ~ M^ 1 ^). This leads to 
an estimate < 2 >~ T M(T), which also may be quite large ||38|| . Here again one should 
be very careful to use the proper coefficients in the estimate. One needs to understand also 
why this estimate could be relevant. 



In order to make the arguments of Ref. [28. 29 more quantitative and to outline the 



domain of their validity, it is helpful to review the stochastic approach to tunneling (see 



39] and references therein). This approach is not as precise as the euclidean approach (in 
theories where the euclidean approach is applicable). However, it is much simpler and more 
intuitive, and it may help us to look from a different point of view on the results we obtained 



in the previous section and on the approach suggested in [28 



The main idea of the stochastic approach can be illustrated by an example of tunneling 
with bubble formation from the point = in the theory (|4lD with the effective potential 

V(<f>, T) = ^f±<f> 2 - £T0 3 + ^0 4 . (48) 

For simplicity, we will study here the limiting case A D — > 0. 

At the moment of its formation, the bubble wall does not move. In the limit of small 
bubble velocity, the equation of motion of the field at finite temperature is simply, 

= d 2 (j)/dr 2 + {2/r)d(j)/dr - V'(<f>) . (49) 

The bubble starts growing if > 0, which requires that 

\d 2 (j)/dr 2 + (2/r)d(j)/dr\ < -V\<f)) . (50) 

A bubble of a classical field is formed only if it contains a sufficiently big field 0. It should 
be over the barrier, so that dV/d<fi < 0, and the effective potential there should be negative 
since otherwise formation of a bubble will be energetically unfavorable. The last condition 
means that the field inside the critical bubble should be somewhat larger than 0*, where 
V(<p*,T) = V(0,T). In the theory (g8|) with < o -> 0, one has 0* = M 2 /2ET. As a simplest 
(but educated) guess, let us take ~ 20* = M 2 /ET. Another important condition is that 
the size of the bubble should be sufficiently large. If the size of the bubble is too small, the 
gradient terms are bigger than the term |V'(0)|, and the field inside the bubble does not 
grow. Typically, the second term in (]50[) somewhat compensates the first one. To make a 
very rough estimate, one may write the condition (]511|) in the form 



\r~ 2 ~ \k 2 < \k 2 max ~ 0- 1 |H0)l ~ 2M 2 . (51) 

Let us estimate the probability of an event in which thermal fluctuations with T ^> M 
build up a configuration of the field satisfying this condition. The dispersion of thermal 
fluctuations of the field with k < k max is given by 



r <^ < ~~.i^max 



k 2 dk 



2 * 2 Jo VW+W fexp ^±f^L 



T r kmax k 2 dk 

(52) 



2tt 2 Jo k 2 + M 2 
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Note that the main contribution to the integral is given by k 2 ~ k 2 max ~ 4M 2 . This means 
that one can get a reasonably good estimate of < 2 >k<k max by omitting M 2 in the integrand. 
This also means that this estimate will be good enough even though the effective mass of 
the scalar field M 2 (0) = V"(4>) changes between = and 0. The result we get is 



T_ [ k ™* _ Tk max C 2 TM 
2Wo 2tt 2 tt 2 



< 2 >*<*_- £j /. dk=±^L = ^ . (53) 



Here C = 0(1) is a coefficient reflecting the uncertainty in the determination of k max and 
estimating the integral. 

Thus, we have a rough estimate of the dispersion of perturbations which may sum up to 
produce a field which satisfies the condition fl51]). We can use it to evaluate the probability 
that these fluctuations build up a bubble of the field of a radius r > k^ ax . This can be 
done with the help of the Gaussian distributionP] 

/ x (f> 2 M 3 ^ 2 , , 4.92M 3 

P{4>) ~ ex P(- 2 <02 >fc<fc _ ) = ex P(-2^PT^ } ~ exp (-^W • ( 54 ) 

Note that the factor in the exponent in ([54]) to within a factor of C 2 = 1.02 coincides with the 
exact result for the tunneling probability in this theory obtained by the euclidean approach 
(see eq. §§)): 

, 4.85M 3 , , , 

P ~ ex P(— &TT) ■ ( 55 ) 

Taking into account the very rough method we used to calculate the dispersion of the per- 
turbations responsible for tunneling, the coincidence is rather impressive. 

As was shown in ||39|| , most of the results concerning tunneling at zero temperature, at a 
finite temperature and even in the inflationary universe, which were obtained by euclidean 
methods, can easily be reproduced (with an accuracy of the coefficient C 2 = 0(1) in the 
exponent) by this simple method. 

Now let us return to the issue of subcritical bubbles. As we have seen, dispersion of 
the long- wave perturbations of the scalar field, < 2 > fc<fcmaa .~ fc ^ 2 T , is quite relevant to 
the theory of tunneling. Its calculation provides a simple and intuitive way to get the same 
results as we obtained earlier by the euclidean approach |39]|. To get a good estimate of the 



probability of formation of a critical bubble in our simple model one should calculate this 
dispersion for k max ~ 2M(T), which gives < 2 >k<k max = TM/tt 2 . Note, that this estimate 
is much smaller than the naive estimate < 2 >~ TM. 

The crucial test of our basic assumptions is a comparison of this dispersion and the value 
of the field at the moment T — T\, when the minimum at = 0i 7^ first appears. Using 

5 The probability distribution is approximately Gaussian even though the effective potential is not purely 
quadratic. The reason is that we were able to neglect the curvature of the effective potential m 2 — V" while 
calculating < 2 >k<k mal - 
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eqs. (I]), (|TTD , one can easily check that the mass of the scalar field at T = T\, — is given 

by 

3ET . . 

m = — =— . (56) 

2y/\ T 



This yields 



For the Higgs boson with m# ~ 60 GeV one obtains 



< 2 > ~ fa « 0! 2_ . (57) 



<h 



< ^ > k<k ~ T • (58) 

Thus, even with account taken of the factor 2/3 in the expression for E, the dispersion of 

long-wave fluctuations of the scalar field is much smaller than the distance between the two 

minima. Therefore, the field on a scale equal to its correlation length ~ M _1 is not equally 

distributed between the two minima of the effective potential. It just fluctuates with a very 

small amplitude near the point = 0. The fraction of the volume of the universe filled by 

the field 0i due to these fluctuations (i.e. due to subcritical bubbles) for m H ~ 60 GeV is 

negligible, 

I <\? \ ( 3En 2 \ _ ia , . 

PW ~ exp (-- < ^ >tt J ~ exp \- w ) ~ e- . (59) 

Since we already successfully applied this method for investigation of tunneling, we expect 
that this estimate is also reliable. The answer remains rather small even for m# ~ 100 GeV, 
when the phase transition is very weakly first order. 

Moreover, even these long-wave fluctuations do not lead to formation of stable domains 
of space filled with the field 0^0, until the temperature is below T c and critical bubbles 
appear. One expects a typical subcritical bubble to collapse in a time r ~ ^mLi this is 
about thirteen orders of magnitude smaller than the total duration of the phase transition, 
At ~ lO -2 // -1 ~ 10~ 4 M p T~ 2 . We do not see any mechanism which might increase r by 
such a large factor. 

Despite all these comments, we think that subcritical bubbles deserve further investi- 
gation. They may lead to interesting effects during phase transitions in GUTs, since the 
difference between T _1 and the duration of the GUT phase transitions is not as great as in 
the electroweak case. They may play an important role in the description of the electroweak 
phase transition as well, in models where the phase transition occurs during a time not much 
longer than T _1 . This may prove to be the case for very weakly first order phase transitions 
with 10 3 GeV ^> niu > 10 2 GeV, when the distance between the two minima of V((j),T) 
at T ~ T\ is smaller than the dispersion y/< 2 >k<k max ~ VTM/tt. 
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7 Conclusions 

One of the main consequences of our work |lj is that it is very difficult to generate baryon 
asymmetry in the standard model without expanding its Higgs sector. One the other hand, 
now we better understand what is necessary for the electroweak baryogenesis to work and 
how to calculate relevant quantities. Hopefully, this will help us to find a realistic theory of 
elementary particles where electroweak baryogenesis is possible. 
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